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11.4 SIGNAL ENHANCEMENT BY TIME-FREQUENCY PEAK 
FILTERING 0 

11.4.1 Signal Enhancement and Filtering 

Time-frequency peak filtering (TFPF) may be regarded as an unconventional alter- 
native to the filtering methods described earlier in this chapter [1] . The signals con- 
sidered are assumed to be sums of arbitrary numbers of band-limited non-stationary 
components in additive noise. For high SNR situations many signal processing al- 
gorithms work well but most perform poorly when SNR decreases below a given 
threshold [2]. In this case, signal enhancement algorithms are required to improve 
the SNR by reducing the distorting effects of noise. To this effect, both adaptive and 
fixed methods have been developed in the case of non-stationary signals in noise. 
Adaptive techniques are generally superior in performance to fixed methods, but 
they perform poorly in certain conditions, such as filtering of a non-stationary sig- 
nal whose spectral content changes rapidly with time. For example, filters designed 
using a least-mean- square (LMS) approach may not adapt quickly enough to track 
the rapidly changing signal due to the delayed convergence of the algorithm. Fur- 
ther, adaptive methods require that the structure of the filter (such as the number 
of the taps) and an estimate of SNR be imposed for optimal performance. This is 
often not possible as assuming a model may lead to suboptimal results and even 
to erroneous conclusions about the signal. This suggests the need for a more gen- 
eral filtering method when the SNR is low and the underlying signal statistics vary 
rapidly with time. 

The TFPF method is based on encoding the noisy signal as the IF of a unit 
amplitude frequency modulated (FM) analytic signal. The instantaneous frequency 
(IF) of the analytic signal is then estimated using standard time- frequency peak 
detection methods [2] to obtain an estimate of the underlying deterministic signal. 
For some signals, TFPF using a windowed WVD results in a significant enhancement 
of signals for SNR as low as —9 dB. 

11.4.2 Time-Frequency Peak Filtering 
11.4.2.1 Background and Definitions 

Let us consider signals expressed as follows: 

p 

s(t) = x(t) + n(t) = ^2 x k{t) + n(t ) (11.4.1) 

k = l 

where n(t) is an additive white Gaussian noise (WGN) and Xk(t) are band-limited 
non-stationary deterministic components that may have overlapping frequency spec- 
tra. It is desired to recover the signal x(t ) given the observation of s(t). 
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The IF of an analytic signal, z(t) 
reference [2]: 



fz(t) 



a(i)ei 2 ' K ^ t \ is defined in Chapter 1 and 



1 d(f>{t) 
2 tt dt 



(11.4.2) 



where <j>(£) is the instantaneous phase and a(t) is the instantaneous amplitude of 
the analytic signal z{t) which can be expressed as [3]: 



z(t) = a(t)e j2n oof-W^ 



(11.4.3) 



Among the existing techniques for IF estimation [2], we adopt the method that 
estimates the IF by taking the peak of the signal’s TFD for its simplicity of imple- 
mentation. The WVD is a natural first choice as a TFD for peak filtering given 
that the other quadratic TFDs are simply smoothed versions of the WVD [3] . The 
WVD of the analytic signal z{t) is defined (in Section 2.1.4) as: 

/ OO 

z{t + r/2 )z*(t — t /2)e _:?27r ^ T dr. (11.4.4) 

-OO 

For monocomponent FM signals, the WVD will produce a time-frequency represen- 
tation of the signal exhibiting significant energy concentration around the signal’s 
IF. When the signal’s IF is linear, delta functions will appear at the positions of the 
IF providing a perfect signal IF estimate. The IF estimate is found by maximizing 
the WVD over frequency [2]; that is 

fz(t) = argmax[IF z (t, /)] (11.4.5) 

/ 

The IF estimate based on the peak of the WVD is unbiased and has variance 
approaching the Cramer- Rao lower bound for signals with linear IF laws in additive 
white zero-mean noise with moderate to high SNR [2]. However, as the order of the 
polynomial IF increases, the delta functions will be replaced by less peaky functions. 
The peak of these functions will lie away from the true IF, resulting in IF estimates 
which are biased. To remedy this, a windowed WVD is used, such that the signal 
IF is as close to linear as possible across the entire window length. 

11.4.2.2 Basic Principle 

Time-frequency peak filtering consists of a two step procedure whereby the signal 
to be filtered is first encoded as the IF of a unit amplitude FM modulated analytic 
signal. Then, the IF is estimated by taking the peak of a time-frequency distribution 
(TFD) to recover the filtered signal. This may be summarized as follows: 

Step 1. Encode the noisy signal s(t) via FM modulation as: z s (t) = e j27r/i ^o s( A) dA 
where /i is a scaling parameter analogous to the FM modulation index. 
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Step 2. Estimate the peak of the WVD of the analytic signal z s (t): 

x(t) = f z (t) = argma x[W Ze (t, f)\//j 
f 

11.4.2.3 Properties 

The properties of the encoding and IF estimation steps of TFPF are derived for the 
case of WVD only. The use of other TFDs will lead to slightly different properties. 

Property 1: The encoding step converts the additive noise n(t) to multiplicative 
noise z n (t) that modulates the signal component z x (t); that is 

z(t) = e j2n »ti s(A)dA = z x (t)z n (t) (11.4.6) 

where the encoded noise and deterministic signal components are given by 

z x (t) and z n (t)=e> 2,rit ti nWdX (11.4.7) 

Proof: Equation 11.4.6 is obtained by a direct substitution of 11.4.1 into 11.4.3. 

Property 2: The Wigner-Ville Spectrum of the signal z s (t) is given by : 

WVS 2s (t, f ) = E[W Zs (: t , /)] = WVS Zn (f , f ) * W 2:c (t, f) (11.4.8) 

where E[-] is the expectation operator, WVS Zn (t,f) = E[W Zn (t, /)], and * repre- 
sents the convolution operation in the frequency domain. 

Proof: This property follows from property 1 and the direct application of the 
expectation operator to the WVD of the encoded signal z s (t). 

This latter is given by: 

E[W Ze (t, /)] = fZ R Zn (t, t)K Zx (t, r)e~^dr (11.4.9) 

where the time-dependent autocorrelation function of z n (t) is 

R Zn (t,T)=E[z n (t + T/2)z*{t-T/2)]=E[e j2 ^ J i-^ nWdX ] (11.4.10) 

and the time-dependent bilinear product function of z x (t) is 

K Zx (t,r) = z x (t + t / 2)z*(t - r/2) = (11.4.11) 

Equation (11.4.8) is then obtained by using the fact that the Fourier transform 
of a product in time is equivalent to the convolution in frequency. 
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Equation (11.4.8) shows that additive noise smears the encoded signal WVD, 
through convolution. Therefore, the bias of TFPF is dependent on 
the encoded signal z x (t) as well as the shape of the encoded noise spectrum 
WVS Zn (t,/). By restricting the shape of the encoded noise spectrum, a class 
of noise can be defined which does not introduce bias to the IF estimation. An 
example from this class is the WGN as will be seen next. 

Property 3: The time dependent autocorrelation function of the encoded noise, 
R Zn (t, r; /i), is equal to the characteristic function of q(t,r) = 2n f*_ ^t/ 2 n iyd\-, 
that is 

R Zn (t,r,y) = E[e? 2 ^^] = $ q (t,T,y) (11.4.12) 

where 4> g (t, r, y) is the characteristic function of q(t,r ) defined by [4]: 

*,(t, r, y) = E[e^o^} = exp ^ (11 .4. 13) 

and k q i(t, r) is the ith cumulant of q(t, r). 

Proof: Equation 11.4.12 is obtained by forming the autocorrelation function of 
z n (t) in equation 11.4.7 and using the above definition of q(t,r). 

11.4.3 Accurate TFPF 

Equation (11.4.8) suggests that in general, a bias in IF estimation is introduced 
by the time-frequency distribution of z x (t ) (deterministic bias) and/or the noise 
(stochastic bias). In the case where the encoded signal s(t) is composed of a deter- 
ministic signal x(t) that is linear in time and embedded in stationary WGN n(t), 
TFPF gives an unbiased estimate of the signal x(t). 

Proof: Consider the signal s(t), given in (11.4.1), to be filtered using TFPF. For 
the case where n(t) is stationary WGN, the ith cumulant of n(t) is such that 
k n i = 0 for i > 3 and q(t, t) is Gaussian with k q i = 0 for i > 3. Furthermore, 
if the noise is a zero-mean independent process, i.e. R n (r) = fc n 2 ^(r), then [5, 
page 369] 

k q i(t, r) = 0 and k q 2 (t, r) = 47r 2 |r|fc„2 (11.4.14) 

The characteristic function given in (11.4.13) becomes 

$,(t, r, y) = e - 2 ’ I ’ 2 ^ 2 l T l fe -2 (11.4.15) 

Taking the Fourier transform of this expression gives 

4n 2 k n2 y 2 

(2Tr 2 k n2 y 2 ) 2 + ( 27 t /) 2 



E[W n (f,t)] 



(11.4.16) 
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This shows that the encoded noise spectrum is low-pass with a maximum at 
the frequency 0 Hz. Hence, WGN will not introduce any bias to the estimate 
of the IF. By replacing this last expression in (11.4.8) we obtain: 



W x .(t,f) = W x ,(t,f)* 



47T 2 hr 



2P 



f (27T 2 fc n2 M 2 ) 2 + (2 7T/) 2 



(11.4.17) 



This expression shows that the bias in the IF using the peak of W Zs (/, t) could 
only come from W Zx (t, /). For the case where the signal x(t) is linear in time; 
that is x(t ) = at + C, where a and C are constants, equation 11.4.17 becomes 



W z .{t,f) = S(f-x(t ))* 



47T 2 /c„2/r 2 



/ (27 r 2 fc n2 M 2 ) 2 
4:TT 2 k n2 p 2 



(2 7 T /) 2 



(27 T 2 k n 2 p 2 ) 2 + (2nf - 2irx(t)) 2 



(11.4.18) 

(11.4.19) 



The delta function ensures that the peak of this function occurs at x(t). There- 
fore if the signal x(t) is linear in time and embedded in stationary WGN, TFPF 
gives an unbiased estimate. 

Equation 11.4.17 shows that in the general case where the signal x(t) is a non- 
linear function of time, the WVD-based TFPF is biased, requiring an appropriate 
windowing of the data. The window is chosen such that the signal within this 
window behaves almost linearly [6] (see Section 11.4.4). In the special case where 
the signal x(t) is a finite-order polynomial in time, the deterministic bias can be 
completely eliminated if the WVD is replaced by the polynomial WVD (PWVD) 
of an appropriate order since the PWVD exhibits delta functions along the IF law 
for polynomial FM signals [7]. 



11.4.4 Discrete-Time Algorithm for TFPF 

The implementation of TFPF using the windowed WVD requires both signal scaling 
before encoding to prevent aliasing, and the selection of the window length for 
reduced bias. These two aspects are discussed next. 



11.4.4.1 Signal Scaling 

FM modulation of un-scaled discrete time signals can lead to aliasing which produces 
discontinuities in the estimated IF at the frequency boundaries of the time- frequency 
plane. This is avoided by amplitude scaling of the noisy signal before frequency 
encoding. Without loss of generality and unless otherwise specified we assume that 
the signal s(t) is sampled at a normalized sampling frequency of 1 Hz. 

The scaled signal, s c (m ), is obtained by using the following transformation. 

s(m) — min[s(m)] 
max[s(m)] — min[s(m)] 



s c (m) = 5[s(m)] = (a — b) 



(11.4.20) 
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where £[.] is the scaling operator and the parameters a and b , which satisfy the 
constraint .5 > a = max[s c (m)] > b = min[s c (m)] > 0, are chosen to provide 
suitable frequency limits on the encoded signal. The operators max[.] and min[.] 
are the maximizing and minimizing functions respectively. The estimate of the 
desired signal, x(m), is recovered by an inverse scaling operation; that is 



x(m) = S 1 [x c (m)\ 



(x c (m) — 6)(max[s(m)] — min[s(ro)]) 
a — b 



+ min[s(m)] 



(11.4.21) 



where x c (m) is the scaled signal obtained using TFPF on s c (m). 



11.4.4.2 Reduced-Bias Window Length Selection 

The bias-variance tradeoff is a key in the practical implementation of TFPF with 
the windowed WVD. Bias reduction requires a small window length to minimize 
the non-optimal nature of the WVD for higher than quadratic phase signals. On 
the other hand, variance reduction is achieved by increasing the window length 
to provide the local estimate with more information. To reduce the variance of 
the estimate while maintaining bias performance it becomes necessary to increase 
the sampling rate. Thus there is a tradeoff between estimator bias and sampling 
frequency which results in a relationship between window length and bias, for a given 
sampling frequency. Results relating to TFPF window lengths are derived in [6]. 
The basic results for window length t w , sampling frequency f s and maximum value 
of IF, f p are given below. For the case of signal estimation: 



t w < 



0.634 f s 
T fp 



(11.4.22) 



This equation gives the maximum window length as a function of maximum signal 
IF and sampling frequency. In a typical application a specified window length is 
required for a given SNR. The signal is sampled at a high enough rate to satisfy 
(11.4.22). 



11.4.4.3 The Iterative TFPF Algorithm 

In the first application, TFPF may not remove as much additive noise as desired. 
If this situation occurs, reapplication of the procedure to the filtered signal is rec- 
ommended. This leads to a 3 step iterative process: 

1. Scale and encode noisy signal. 

2. Apply TFPF to yield a signal estimate x{t). 

3. If x(t) contains substantial noise, go to step 1; else terminate the process. 

Simulations demonstrate the convergence of the repeated scheme to a stable esti- 
mate x(t) [1]. 
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(a) Original synthetic signal (d) Original EEG signal 




(b) Noisy synthetic signal (e) Noisy EEG signal. 




(c) Enhanced synthetic signal (f) Enhanced EEG signal 

Fig. 11.4.1: B-distributions of a synthetic signal (left column) and an EEG signal (right column), 
showing the original signal (top row), and the noisy signal before enhancement (middle row) and after 
enhancement (bottom row). 
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11.4.5 Examples and Results 

Example 1 (A Multicomponent signal in WGN) Let us consider a multicom- 
ponent signal, x(m), expressed by: 



x(m) 



' 0.85 sin(0.055m + 3.75 x 10~ 4 sin(0.000625m)m) + sin(0.035m) 

+ sin(0.020m), 0 < m < L/2 

< 0.85 sin(0.055m + 3.75 x 10~ 4 sin(0.000625m)m) 

+2(1 - m/L) sin(0.035m + 2.75 x 10- 11 (m - L/ 2) 3 ) 

+2(1 - m/L ) sin(0.020m - 2.75 x Kr n (m - L/ 2) 3 ), L/2<m<L 



where the data length L is taken as 32768 data points. For a time-frequency illus- 
tration of this signal, the B-distribution (BD) with smoothing parameter /3 = 0.01 
is shown in Fig. 11.4.5(a). 1 White Gaussian noise was added to the above signal 
giving an SNR of — 9dB; the BD of the noisy signal is shown in Fig. 11.4.5(b). The 
windowed WVD peak filter was then implemented to recover x[m ) from the noisy 
signal. A window length of 15 data points was chosen to satisfy the window length 
constraints given in (11.4.22). Fig. 11.4.5(c) shows the clean recovery of the signal 
after three TFPF iterations. Note that the WVD is used as the vehicle for signal 
recovery while the B-distribution is used only for presentation of the results. 

Example 2 (Newborn EEG data in WGN): Fig. 11.4.5(d) shows a time- 

frequency representation of a real newborn EEG signal using the B-distribution 
with /3 = 0.01. WGN is then added to the signal at SNR = —9 dB. The noisy 
signal in Fig. 11.4.5(e) shows that the time-frequency patterns of the EEG signal 
are not clearly visible. Using a window length of 20 data points, four iterations of 
the TFPF were used to recover a cleaner signal. The filtered signal in Fig. 11.4.5(f) 
demonstrates the efficiency of TFPF. 



11.4.6 Summary and Conclusions 

TFPF is a tool for signal enhancement, applicable to a large class of signals if 
the windowed WVD TFPF is used for reduced bias. This class includes those 
signals which may be represented as a sum of band-limited non-stationary processes 
in additive WGN. Testing on simulated and real data indicates that the method 
significantly enhances signals of this class by filtering out most of the additive noise. 
Further details of the time- frequency peak filtering method are provided in [8] . 
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